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1. Introduction 

In this topical review we briefly discuss the status of astronomical black hole 
observations and focus on phenomenology of the mini-black- hole production and decay 
in high-energy collisions in models with low-scale gravity. We point out exciting ways 
of studying quantum gravity and searching for new physics using large samples of 
black holes that may be accessible at future colliders and discuss the potential of 
the existing and future cosmic ray detectors for searches for black hole production in 
ultra-high-energy cosmic rays. 

In what follows we will use "natural" units: 

7i = c = Aib = 1, 

where h is the reduced Planck's constant, c is the speed of light, and fcs is Boltzmann's 
constant. This choice allows us to measure both energy and temperature in the 
same units, TeV (= 10^^ eV), while distance and time are measured in TeV~^. 
The relationship between natural and SI units is as follows: 1 TeV= 1.16 x 10^^ K, 
1 TeV"^ = 6.58 x 10"^^ s = 1.97 x IQ-^^ m. 

2. The hierarchy problem 

One of the most pressing problems in modern particle physics is the hierarchy problem, 
which can be expressed by a single, deceptively simple question: why is gravity (at least 
as we observe it) some 10^® orders of magnitude weaker than other forces of Nature? 
This mysterious fact is responsible for a humongous hierarchy of energy scales: from 
the scale of electroweak symmetry breaking (Mew ~ 1 TeV) to the energies at which 
gravity is expected to become as strong as other three forces (strong, electromagnetic, 
and weak), known as the Planck scale, or Mpi = l/y'Gjv ~ 1-22 x 10^^ TeV, where 
Gjv is Newton's coupling constant. 

If it was not for the hierarchy problem, our current understanding of fundamental 
particles and forces, consolidated in the standard model (SM) of particle physics, would 
have been nearly perfect. After all, the standard model has an impressive calculational 
power. Its predictions have been tested with a per mil or better accuracy, with hardly 
any significant deviations found so far. While certain phenomena, such as neutrino 
masses or existence of dark matter, are not explained in the standard model, they 
still may be accommodated with its minimum expansion, just like CP- violation has 
been earlier added to the SM framework via the quark mixing matrix. While it is true 
that several phenomena, including CP-violation or the hierarchy of fermion masses, 
are not really explained within the standard model, nevertheless they are readily 
accommodated within the SM by introducing additional parameters, e.g. fermion 
couplings to the Higgs field. Even the only remaining unobserved particle predicted 
in the standard model a fundamental Higgs boson responsible for the weak boson 
masses, while leaving photon massless - may soon be found at the Large Hadron 
Collider (LHC, currently under construction at CERN, near Geneva, Switzerland) 
or even at the Tevatron (Fermilab, near Chicago, USA). Consequently, there are no 
fundamental experimental reasons to believe that the standard model cannot work 
to the highest energies, such as Planck scale. Nevertheless, there is a pretty good 
theoretical reason to believe that this is not the case: the hierarchy problem. 

The reason that the hierarchy problem is so annoying stems from the fact that in 
our everyday experience large hierarchies tend to collapse, unless they are supported 
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by some intermediate layers or by precisely tuned initial conditions (the so-called 
"fine tuning"). For example, while it is mechanically possible to balance a pen on 
its point, the amount of fine tuning required for such a delicate balance is simply 
too demanding to be of practical use. Indeed, walking in a room and seeing a pen 
standing vertically on its point on a table would be too odd not to suspect some kind 
of a hidden support. Similarly, making the standard model work for all energies up to 
the Planck scale would require tremendous amount of fine tuning of its parameters - 
to a precision of ^ (Mpi/Mew)^ ^ 10"'^^! A natural reaction to this observation is 
to conclude that the standard model needs some "support," which should come either 
in a form of new physics at intermediate energy scales that adds extra layers to the 
hierarchy, or from certain new symmetries that would guarantee necessary amount of 
fine tuning. 

However, it is worth pointing out that large amounts of fine tuning, while 
unnatural, are not prohibited by any fundamental principle. The fact that large 
hierarchies tend to collapse is merely empirical. While many large hierarchies we know 
have indeed collapsed, some other are surprisingly stable - examples of both kinds can 
be found in mechanical, social, and political systems alike. Certain examples of fine 
tuning have been observed in Nature: from the fact that the apparent angular size of 
the moon is the same as the angular size of the Sun within 2.5% - a mere coincidence 
to which we owe such a spectacular sight as a solar eclipse to a somewhat less abused 
example of the infamous Florida 2000 presidential election recount with the ratio of 
Republican to Democratic votes equal to 1.000061, i.e. fine-tuned to unity with the 
precision of 0.006%! 

Despite these observations, the majority of physicists have been working very 
hard to find a more rational explanation of the hierarchy of forces in the standard 
model since its formulation in the late sixties. A number of viable solutions have been 
proposed as a result of this work: from supersymmetry (which protects the hierarchy 
due to nearly exact intrinsic cancelations of the effects caused by the standard model 
particles and by their superpartncrs obeying different spin statistics) to models with 
strong dynamics (which introduce an intermediate energy scale via new, QCD-like 
force, and often result in a composite Higgs boson). 

Staying somewhat aside and quite debatable and controversial is an explanation 
of the standard model hierarchy via the anthropic principle: our universe is so fine 
tuned, because it happens to be one of a very few possible universes that could support 
intelligent life capable of raising this very question. Sparkled by the recent findings 
that the number of possible string theory vacua may be tremendously large: say, 10^°° 
or may be even many orders of magnitude more, the anthropists claim that this vast 
landscape can give rise to enormous variety of possible universes, thus making it very 
likely the existence of significantly fine-tuned ones, which could support formation of 
heavy elements, stars, galaxies, and other precursors of intelligent life. For discussion 
of the recent controversy around the anthropic principle, sec [1]. 

While each of the aforementioned solutions to the hierarchy problem easily 
deserves a separate review, here we will focus on yet another, more recently proposed 
remedy that does not involve new particles or symmetries, but instead uses the 
geometry of space itself. 
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3. A brief history of space 

The idea that our space may contain more than the three famiUar dimensions has been 
one of the most popular recurring themes in the work of numerous philosophers, artists, 
and writers since ancient times. From the very concept of Heaven and Hell to shadow 
or mirror worlds and parallel universes, the believe that the space around us may be 
a bigger place than it is commonly thought became one of the most mesmerizing and 
puzzling ideas, which inspired many great minds of the past centuries. However, it 
took quite a while before this concept has been truly embraced by scientists. 

In fact, in the beginning of the nineteenth century, following the seminal work 
by Gauss, there have been numerous publications that claimed that additional 
dimensions in space would contradict the inverse square law obeyed by gravitational 
and electromagnetic forces. Consequently the very concept of extra dimensions was 
quickly abandoned by mathematicians and physicists and remained the realm of art 
and philosophy. 

Perhaps the first person to consider the possibility of extra dimensions in space 

rigorously was Bernhard Riemann [2]. What started as an abstract mathematical 
idea of a curved Riemannian space, soon became the foundation of the most profound 
physics theory of the last century, if not of the entire history of physics: Albert 
Einstein's general relativity [3]. While Einstein's theory was formulated in the three- 
plus-one space-time dimensions, it soon became apparent that the theory cannot be 
self-consistent up to the highest energies in its original form. 

In the 1920s, Theodor Kaluza and Oskar Klein [4] suggested that a unification of 
electromagnetism and general relativity is possible if the fifth, spatial dimension of a 
finite size is added to the four-dimensional space-time. While this attempt has not led 
to a satisfactory and self consistent unification of gravity and electromagnetism, the 
idea of finite (or "compactified" ) extra spatial dimensions has been firmly established 
by Kaluza and Klein and eventually led to their broad use in string theory. Rapid 
progress in string theory in the 1970s helped the original idea of half-a-century earlier 
to regain its appeal. It was realized that extra six or seven spatial dimensions are 
required for the most economical and symmetric formulation of string theory. In 
particular, string theory requires extra dimensions to establish its deep connection 
with the supersymmetry, which also leads to the unification of gauge forces. While 
the size of these compact dimensions is not fixed in string theory, it is natural to 
expect them to have the radii similar to the inverse of the grand unification energy 
scale, or of the order of 10^'^^ m. Unfortunately, no experimental means exist to probe 
such short distances, so extra dimensions of string theory will likely remain untested 
even if this theory turns out to be the correct description of quantum gravity. 

4. Large extra spatial dimensions 

The situation with seemingly untestable extra spatial dimensions changed dramatically 
in the late 1990s, when Arkani-Hamed, Dimopoulos, and Dvali (ADD) suggested a 
new paradigm [5] in which several (n) of compactified extra dimensions could be as 
large as ~ 1 mm! These large extra dimensions have been introduced to solve the 
hierarchy problem of the standard model by dramatically lowering the Planck scale 
from its apparent value of Mpi ^ 10^^ TeV to ^ 1 TeV. (We further refer to this 
low fundamental Planck scale as Mz).) In the ADD model, the apparent Planck scale 
Mpi only reflects the strength of gravity from the point of view of a three-dimensional 
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(3D) observer. It is, in essence, just a virtual "image" of the fundamental, (3 + n)- 
dimcnsional Planck scale Mjj, caused by an incorrect interpolation of the gravitational 
coupling (measured only at low energies and large distances) to a completely different 
regime of high energies and short distances. 

In order for such "large" extra dimensions not to violate any constraints from 
atomic physics and other experimental data, all other forces except for gravity must 
not be allowed to propagate in extra dimensions. Conveniently, modern quantum field 
theory allows to confine spin 0, 1/2, and 1 particles to a subset of a multidimensional 
space, or a "brane." At the same time, the graviton, being a spin-2 particle, is not 
confined to the brane and thus permeates the entire ("bulk") space. The solution 
to the hierarchy problem then becomes straightforward: if the Planck scale in the 
multidimensional space is of the order of the electroweak scale, there is no large 
hierarchy and no fine-tuning problem to deal with, as non-trivial physics really "ends" 
at the energies ^ 1 TcV. Gravity appears weak to a 3D-obscrvcr merely because it 
spans the entire bulk space, so the gravitons spend very little time in the vicinity 
of our brane. Consequently, the strength of gravity in 3D is literally diluted by an 
enormous volume of extra space. 

The compactification radius of large extra dimensions (i?c) is fixed by their 
number and the value of the fundamental Planck scale in the (4 -|- n)-dimensional 
space-time. By applying Gauss's law, one finds [5] 

M|i = 87rM^+2i?J^. (1) 

Here for simplicity we assumed that all n extra dimensions have the same size Rc and 
are compactified on a torus. 

If one requires Md ^ I TeV and a single extra dimension, its size has to be as 
large as the radius of the solar system; however, already for two extra dimensions their 
size is just ~ 1 mm; for n = 3 it is ~ 1 nm, i.e. similar to the size of an atom; and for 
larger n the radius further decreases to a subatomic size and reaches 1 fm (the size 
of a proton) for seven extra dimensions. 

A direct consequence of such compact extra dimensions in which gravity is allowed 
to propagate, is a modification of Newton's law at the distances comparable with Rc'. 
the gravitational potential would fall off as l/r"+^ for r < Rc- This immediately rules 
out the possibility of a single extra dimension, as the very existence of our solar system 
requires the potential to fall off as l/r at the distances comparable with the size of 
planetary orbits. However, as of 1998 any higher number of large extra dimensions 
have not been ruled out by gravitational measurements, as Newton's law has not been 
tested to the distances smaller than about 1 mm. As amazing as it sounds, large extra 
dimensions with the size as macroscopic as ^ 1 mm were perfectly consistent with the 
host of experimental measurements at the time when the original idea has appeared! 

Note that the solution to the hierarchy problem via large extra dimensions is not 
quite rigorous. In a sense, the hierarchy of energy scales is traded for the hierarchy 
of distance scales from the "natural" electroweak symmetry breaking range of 1 
TeV^^ ^ 10^^® m to the much larger size of compact extra dimensions. Nevertheless, 
this hierarchy could be significantly smaller than the hierarchy of energy scales or 
perhaps even be stabilized by some topological means. Strictly speaking, the idea of 
large extra dimensions is really a paradigm, which can be used to build more or less 
realistic models, rather than a model by itself. 

Since large extra dimensions as a possible solution for the hierarchy problem 
were suggested, numerous attempts to either find them or to rule out the model 
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have been carried out. They include measurements of gravity at sub-millimeter 
distances [6], studies of various astrophysical and cosmological implications [7], and 
numerous collider searches for virtual and real graviton effects [8]. For a detailed 
review of the existing constraints and sensitivity of future experiments, the reader is 
referred to [8, 9, 10]. The host of experimental measurements to date have largely 
disfavoured the case of two large extra dimensions. However for three or more extra 
dimensions, the lower limits on the fundamental Planck scale are only 1 TeV, i.e. 
most of its allowed range has not been probed experimentally yet. 

5. The Randall- Sundrum model 

A significantly different and perhaps more rigorous solution to the hierarchy problem 
is offered in the Randall-Sundrum (RS) model [11] with non-factorizable geometry 
and a single compact extra dimension. This is achieved by placing two 3-dimensional 
branes with equal and opposite tensions at the fixed points of the S1/Z2 orbifold in 
the five-dimensional anti-deSitter space-time The metric of the AdS^ space is 

given by rfs^ = cxp(—2kRc\Lp\)r]^,ydx'^dx'^ — Blf.d'.f^ , where < < tt is the coordinate 
along the compact dimension of radius Rc, k is the curvature of the AdS5 space, often 
referred to as the warp factor, are the conventional (3+l)-space-time coordinates, 
and r]'^'^ is the metric tensor of Minkowski space-time. The positive-tension (Planck) 
brane is placed at = 0, while the second, negative-tension standard model brane is 
placed at <f = TV. In this framework, gravity originates on the Planck brane and the 
graviton wave function is exponentially suppressed away from the brane along the fifth 
dimension due to the warp factor. Consequently, the O(Mpi) operators generated on 
the Planck brane yield low-energy effects on the standard model brane with a typical 
scale of = Mp\ exp(— 7rfci?c), where Mpi = A'/pi/a/Stt is the reduced Planck mass. 
The hierarchy between the Planck and electroweak scales is solved if Att is ~ 1 TeV, 
which can be achieved with little fine tuning by requiring kRc ~ 12. It has been 
shown [12] that the size of the extra dimension can be stabilized by the presence of 
a bulk scalar field (the radion). Consequently, the hierarchy problem in the Randall- 
Sundrum model is solved naturally for k ~ 10/ Rc ^ (10~^ — 10~^)Mpi, since Mp\ is 
the only fundamental scale in this model and both k and 1/Rc are of the same order 
as this scale. 

It is convenient to introduce a dimensionless parameter k = k/Mpi, which 
defines the strength of coupling between the graviton and the standard model fields. 
Theoretically preferred value of k is between 0.01 and 0.1. For larger coupling theory 
becomes non-perturbative; if the coupling is too small, an undesirably high amount 
of fine tuning is still required to solve the hierarchy problem. 

In the simplest form of the Randall-Sundrum model only gravitons are allowed 
to propagate in the fifth dimension; in the extensions of the model other particles 
are allowed to "live" in the bulk [13]. In both cases, since gravitons are allowed to 
propagate in the bulk, from the point of view of a 4D-observer on the standard model 
brane, they acquire Kaluza-Klein modes with the masses given by the subsequent 
zeroes Xi of the Bessel function Ji {Ji{xi) — 0): Mi = kxie~'^^^'' = kxiAT^- The mass 
of the first excited mode is then Mi « 3.83A;A^ ~ 1 TeV. The zeroth Kaluza-Klein 
mode of the graviton remains massless and couples to the standard model fields with 
the gravitational strength, 1/Mpi, while the excited graviton modes couple with the 
strength of l/A^r. Coupling k governs both the graviton width and the production 
cross section at colliders (both are proportional to fc^). 
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Figure 1. Limits on the Randall-Sundrum model parameters Mi and k = k/Mpi 
from tlie D0 experiment [14]. The light-shaded area has been excluded in the 
dimuon channel; the medium-shaded area shows the exclusion obtained in the 
dielectron and diphoton channels; the dotted line corresponds to the combination 
of all three channels. The area below the dashed-dotted line is excluded from the 
precision electroweak data (see [13]). The dark shaded area in the lower right- 
hand corner corresponds to A.^- > 10 TeV, which requires undesirably high amount 
of fine tuning. 



Indirect limits on the Randall-Sundrum model parameters come from precision 
electroweak data (dominated by the constraint on the S parameter) [13]. The 
only dedicated search for Randall-Sundrum gravitons at colliders so far has been 
accomplished by the D0 experiment at the Fermilab Tevatron and set direct limits [14] 
on A: as a function of the mass of the first Kaluza-Klein mode of graviton, Mi. Both 
direct and indirect limits are summarized in figure 1. As seen from the figure, for 
masses Mi w 3.83fcA7r > 800 GeV most of theoretically preferred range for k is still 
allowed. 



6. Astronomical black holes 

While very few scientists doubt that black holes exist somewhere in the universe, and 
perhaps are even abundant, all the astronomical black hole candidates discovered so 
far are only found via indirect means. There are several ways astronomers look for 
black holes. For example, one of the first established black hole candidates, Cygnus 
X-1 binary system, was found by observing the orbital periods in superbright binary 
systems. The presence of a black hole as one of the two stars in Cygnus X-1 was inferred 
from the large total power dissipated by the system, large estimated mass of the 



Black holes at future colliders and beyond 



8 



unseen companion star, and from extremely short time-scale of the intensity variations. 
Similar arguments load astronomers to believe that quasars are powered by massive 
black holes. Furthermore, by observing X-ray flares in the active galactic nuclei, it 
is speculated that they are caused by large objects falling inside the nuclei, being 
attracted by supermassive 10^ solar masses) black holes located in their centers. 
The presence of these supermassive black holes is also inferred by the motion of stars 
near the galactic centers. There is even an evidence for such a supermassive black 
hole in the center of our own galaxy, the Milky Way. Clearly, these observations prove 
the presence of massive compact objects in many of the binary systems or galaxies; 
however, proving that the critical density, necessary for a gravitational collapse into 
a black hole, has been achieved in these systems is rather complicated. The recent 
announcement of an observation of a "strange" star, apparently more dense than a 
neutron star [15], given the lack of a compelling cosmological explanation or even 
motivation for such an object, might indicate that there is a problem with the current 
methods of density estimates in extremely remote systems. 

Even the direct dynamical observations only probe the region of space at least two 
orders of magnitude further away from the center of a black hole candidate than the 
estimated radius of its event horizon (often referred to as the Schwarzschild radius, 
Rs)- To date, the closest we have been able to probe our own galactic center is 
w 600i?s [16]. 

Other phenomena predicted for black holes in general relativity, such as frame 
dragging, have been observed as well. However, none of the existing black hole 
candidates have been tagged by several independent means so far. Given the large 
number of objects studied by the astronomers in their quest for black holes, it is 
questionable how unambiguous the single tags are. 

Unfortunately, the most prominent feature of a black hole — its Hawking 
radiation [17] — has not been observed yet and is very unlikely to be ever observed 
by astronomical means. Indeed, even the smallest (and therefore the hottest) 
astronomical black holes with the mass close to the Oppenheimer-Volkov limit [18] 
of ~ 3 solar masses, have Hawking temperatures of only ~ 100 nK, which corresponds 
to the wavelength of Hawking radiation of ~ 100 km, and the total dissipating power 
of puny - 10^28 ^1 

Not only the event horizon of these black holes is nearly as cold as the lowest 
temperature ever achieved in the lab (the 1997 Nobel prize in physics [19] was given 
to Chu, Cohen- Tannoudji, and Phillips who reached the temperature as low as ^ 1 /iK 
via optical cooling), but the wavelength of the radiation dissipated by a black hole is 
far from the visible spectrum and is resemblant of that of an AM radio station. Trying 
to detect such a radio broadcast with a vanishing transmitting power from thousands 
of light-years away is but impossible. Given that the black hole dissipating power 
corresponds to only ~ 100 photons per second emitted by its entire event horizon and 
that the closest known black hole candidate is still over a thousand of light years away 
from us, not a single Hawking radiation photon ever hit our Earth since it has been 
formed! (In fact, one would have to wait 10^^ years to observe a single photon from 
such a black hole to hit the Earth.) Thus, if the astronomical black holes were the 
only ones to exist, the Hawking radiation would be always just a theoretical concept, 

X Note that the event horizon temperature of these black holes is much lower than the temperature 
of the cosmic microwave background radiation, so at the present time the black holes are growing due 
to the accretion of relic radiation much faster than they are evaporating. The evaporation process 
will take over only when the expanding universe cools down below the event horizon temperature. 
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never testable experimentally. 

While Hawking radiation would constitute a definitive proof of the black hole 
nature of a compact object, there are other, indirect means of identifying the existence 
of the event horizon around it. It has been suggested [20] that the lack of Type I X-ray 
bursts in the binary systems identified as black hole candidates, implies the presence 
of the event horizon. The argument is based on the fact that such X-ray bursts 
are frequent in the neutron star binaries, similar in size and magnitude to the black 
hole candidates. Consequently, if the black hole candidates did not have the event 
horizon, one would expect to see a similar X-ray burst activity, which contradicts 
the observations. An analogous argument applies to the X-ray supernovae that are 
believed to contain black holes. These supernovae are much dimmer than similar 
supernovae that are believed to contain neutron stars, which is considered to be an 
evidence for the formation of an event horizon around the black hole candidates. 
Unfortunately, an evidence based on a non-observation of a particular predicted 
phenomenon is inherently much more model-dependent and circumstantial than the 
one based on an observation of a certain effect. Perhaps, a more promising way to 
prove the existence of the event horizon around some of the black hole candidates is to 
compare the accretion disk shapes in various X-ray binaries [21]. Black hole binaries 
are expected to have advection-dominated accretion flows, drastically different from 
thin accretion disks, typical of subcritical binaries. 

Probably, the best evidence for the existence of astronomical black holes would 
come from an observation of gravitational waves created in the collisions of two black 
holes, which LIGO and VIRGO detectors are looking for. However, current sensitivity 
of these interferometers is still short of the expected signal, even in optimistic 
cosmological scenarios. 

This leads us to other places to look for black holes that are much smaller and 
consequently much hotter and easier to detect than their astronomical counterparts. 

7. Mini black holes in large extra dimensions 

As was pointed out several years ago [22, 23, 24], an exciting consequence of low- 
scak; quantum gravity is the possibility of producing black holes (BH) at particle 
accelerators. More recently, this phenomenon has been quantified [25, 26]for the case 
of particle collisions at Large Hadron Collider (LHC), resulting in a mesmerizing 
and unexpected prediction that future colliders would produce mini black holes at 
enormous rates (e.g., ~ 1 Hz at the LHC for Md = 1 TeV), thus qualifying for black- 
hole factories. With the citation index of the original papers [25, 26] well over three 
hundred, the production of mini black holes in the lab became one of the most actively 
studied and rapidly evolving subjects in the phenomenology of models with low-scale 
gravity 

In general relativity, black holes are well understood if their mass A'/bh far exceeds 
the Planck scale. Consequently, in the model with large extra dimensions general 
relativity would give accurate description of black-hole properties when its mass is 
much greater than the fundamental (multidimensional) Planck scale Mo ^ 1 TeV. 
As its mass decreases and approaches Md, the black hole becomes a quantum gravity 
object with unknown and presumably complex properties. 

In this section, following [25], we will ignore this obstacle§ and estimate the 

§ Some of the properties of the "stringy" subplanclfian "precursors" of blaclj holes are discussed in 
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properties of light black holes by simple semiclassical arguments, strictly valid only 
for Mbh ^ Mjj. Wc expect this to be an adequate approximation, since the important 
experimental signatures rely on two simple qualitative properties: the absence of small 
couplings and the "democratic" nature of black hole decays, both of which may survive 
as average properties of the light descendants of black holes. 

As we expect unknown quantum gravity effects to play an increasingly important 
role for the black hole mass approaching the fundamental Planck scale, following the 
prescription of [25], we do not consider black hole masses below the Planck scale. It 
is expected that the black-hole production rapidly turns on, once the relevant energy 
threshold ~ Md is crossed. At lower energies, we expect black-hole production to be 
exponentially suppressed due to the string excitations or other quantum effects. 

We will first focus on the production in particle collisions and subsequent decay 
of small Schwarzschild black holes with the size much less than the compactification 
radius of extra dimensions. In this case, standard Schwarzschild solution found for 
a flat (3 + n)-dimensional metric fully applies. The expression for the Schwarzschild 
radius Rs of such a black hole in (3 + n) spacial dimensions is well known [28]: 



Rs{Mbh) = ^ 



Mbh SF (^±2) 



Md n + 2 



(2) 



Given Mjj ~ 1 TeV and taking into account the fact that black hole masses accessible 
at the next generation of particle colliders and in ultrahigh-energy cosmic ray collisions 
are at most a few TeV, we note that the Schwarzschild radius of such black holes is 
~ l/Mjj, i.e. indeed much smaller than the size of large extra dimensions even 
when their number approaches six or seven (the preferred number of extra dimensions 
expected in string theory). We also note that for Mbh ~ Md the Schwarzschild radius 
does not depend significantly on the number of extra dimensions n. 

Given the current lower constraints on the fundamental Planck scale in the model 
with large extra dimensions of ~ 1 TeV [8], the black holes that we may be able 
to study at colliders and in cosmic rays will be barely transplanckian. Hence, the 
unknown quantum corrections to their classical properties are expected to be large, 
and therefore it is reasonable to focus only on the most robust properties of these mini 
black holes that are expected to be affected the least by unknown quantum gravity 
corrections. Consequently, when discussing production and decay of black holes, we 
do not consider the effects of spin and other black hole quantum numbers, as well as 
grey-body factors. Further in the review, we discuss some of the subsequent attempts 
to take these effects into account using semiclassical approximation. 



7.1. Black-hole production in particle collisions 

Consider two partons with the center-of-mass energy \f§ — Mbh colliding head-on. 
Semiclassical reasoning suggests that if the impact parameter of the collision is less 
than the Schwarzschild radius Rs, corresponding to this energy, a black hole with the 
mass Mbh is formed. Therefore the total cross section of black hole production in 
particle collisions can be estimated from pure geometrical arguments and is of order 
7rJ?|. 

Soon after the original calculations [25, 26] have appeared, it has been 
suggested [29] that the geometrical cross section is in fact exponentially suppressed, 

[27] and later in this review. 
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Figure 2. a) Parton-lcvcl and b) differential production cross section of black 
holes in pp collisions at the LHC. Also shown: c) Hawking temperature and d) 
average particle multiplicity in black-hole decays. Adapted from [25]. 



based on the Gibbons-Hawking action [30] argument. Detailed subsequent 
studies performed in simple string theory models [27], using full general relativity 
calculations [31], or a path integral approach [32] did not confirm this finding and 
proved that the geometrical cross section is modified only by a numeric factor of order 
one. A flaw in the Gibbons-Hawking action argument of [29] was further fotmd in [33]: 
the use of this action implies that the black hole has been already formed, so describing 
the evolution of the two colliding particles before they cross the event horizon and form 
the black hole via Gibbons-Hawking action is not justified. By now there is a broad 
agreement that the production cross section is not significantly suppressed compared 
to a simple geometrical approximation, which we will consequently use through this 
review. 

Using the expression (2) for the Schwarzschild radius [28] , we derive the following 
parton level black hole production cross section [25] : 

2 

n + 1 

(3) 



(7(Mbh) « ttRI = 



M] 



BH 



D 



n + 2 



This cross section, as a function of the fundamental Planck scale and the black hole 
mass, is shown in figure 2a. 
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In order to obtain the production cross section in pp collisions at the LHC, we 
use the parton luminosity approach [25, 26, 34]: 



dcT{pp^BR + X) dL 

— ^ Tin = -m — ^(«^ BH) 

dMBH dMBH 



s=Mi 



where the parton luminosity dL/dM-QB is defined as the sum over all the types of 
initial partons: 



dL _ 2Mbh \p / axa ^ ^^.^ ivi^u ^ 



a,b "^"bh/ 

and fi{xi) arc the parton distribution functions (PDFs). We used the MRSD— ' [35] 
PDFs with the scale taken to be equal to Mbh, which is within the allowed range 
of these PDFs for up to the kinematic limit at the LHC. The dependence of the cross 
section on the choice of PDF is only ^ 10%. The total production cross section for 
Mbh > Md at the LHC, obtained from the above equation, ranges between 15 nb 
and 1 pb for the Planck scale between 1 TeV and 5 TeV, and varies by 10% for n 
between 2 and 7. The differential cross section is shown in figure 2b. 

7.2. Black-hole evaporation 

In general relativity, black-hole evaporation is expected to occur in three distinct 
stages: "balding," spin-down, and Hawking evaporation. During the first stage, the 
black hole loses its multipole momenta and quantum numbers via emission of gauge 
bosons until it reaches the Kerr solution for a spinning black hole; at the second stage 
it gets rid of the residual angular momentum and becomes a Schwarzschild black 
hole; and at the last stage it decays via emission of black-body radiation [17] with a 
characteristic Hawking temperature: 

1 

f Md n + 2 \"^' n+1 n + 1 

^" = '''' [m^WW) J V^ = i^ ^ ^ 

of ~ 1 TeV [25, 26]. The dependence of the Hawking temperature on the fundamental 
Planck scale and the black-hole mass is shown in figure 2c. 

Note that if a certain quantum number (e.g., B — L) is gauged, it will be conserved 
in the process of black hole evaporation. Since the majority of black holes at the LHC 
are produced in quark-quark collisions, one would expect many of them to have the 
baryon number and fractional electric charge. Consequently, the details of black- 
hole evaporation process will allow to determine if these quantum numbers are truly 
conserved. 

In quantum gravity, it is expected that there is a fourth, Planckian stage of 
black hole evaporation, which is reached when the mass of the evaporating black 
hole approaches the Planck scale. The details of the Planckian stage are completely 
unknown, as they are governed by the effects of quantum gravity, which should be 
dominant at such low black hole masses. Some authors speculate that the Planckian 
stage terminates with a formation of a stable or semi-stable black hole remnant with 
the mass ^ M-pi. Others argue that the evaporation proceeds imtil the entire mass 
of the black hole is radiated. The truth is that no predictions about the Planckian 
regime are possible, given our lack of knowledge of quantum gravity. 
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The average multiplicity of particles produced in the process of black-hole 
evaporation is given by: (N) = ^ MsR ^ ^ where E is the energy spectrum of the decay 

products. In order to find (A^), wc note that evaporation is a black-body radiation 
process, with the energy flux per unit of time given by Planck's formula: 

^ ~ (5) 
dx + c ' 

where x = E/Th, and c is a constant, which depends on the quantum statistics of the 
decay products (c = —1 for bosons, +1 for fermions, and for Boltzmann statistics). 

The spectrum of the decay products in the massless particle approximation is 
given by: ^ ~ ^ ~ j^fj • For averaging the multiplicity, we use the average of 
the distribution in the inverse particle energy: 

* ' JO e^+c 

where a is a dimensionless constant that depends on the type of produced particles and 
numerically equals 0.68 for bosons, 0.46 for fermions, and | for Boltzmann statistics. 
Since a mixture of fermions and bosons is produced in the black hole decay, we can 
approximate the average by using Boltzmann statistics, which gives the following 
formula for the average multiphcity: {N) ^ Using expression (4) for Hawking 

temperature, we obtain: 




which corresponds to about half-a-dozen for typical black hole masses accessible at 
the LHC, see figure 2d. 

Naively, one would expect that a large fraction of Hawking radiation is emitted 
in the form of gravitons, escaping in the bulk space. However, as was shown in [24], 
this is not the case, since the wavelength A = ^ corresponding to the Hawking 
temperature is larger than the size of the black hole. Therefore, the black hole acts as 
a point-radiator and consequently emits mostly s-waves. Since the s-wave emission is 
sensitive only to the radial coordinate, bulk radiation per graviton degree of freedom 
is the same as radiation of any standard model degree of freedom on the brane. While 
many angular degrees of freedom are available in the bulk space, the s-wave emission 
cannot take advantage of them, thus suppressing bulk graviton component. Since 
there are many more particles on the brane than in the bulk space, this has the crucial 
consequence that the black hole mainly decays to visible standard model particles. 

Since the gravitational coupling is flavour-blind, a black hole emits all the w 120 
standard model particle and antiparticle degrees of freedom with roughly equal 
probability. Accounting for the colour and spin and ignoring the graviton emission, 
we expect « 75% of particles produced in black hole decays to be quarks and gluons, 
w 10% charged leptons, « 5% neutrinos, and « 5% photons or W/Z bosons, each 
carrying hundreds of GeV of energy. Similarly, if there exist new particles with masses 
< 100 GeV, they would be produced in the decays of black holes with the probability 
similar to that for the standard model species. For example, a sufficiently light Higgs 
boson is expected to be emitted with ^1% probability. This has exciting consequences 
for searches for new physics at the LHC and beyond, as the production cross section for 
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Figure 3. Number of black holes produced at the LHC in the electron or photon 
decay channels (assuming 15% efficiency for an electron or a photon tag), with 
100 fb~^of integrated luminosity, as a function of the black hole mass. The shaded 
regions correspond to the variation in the number of events for n between 2 and 
7. The dashed line shows total standard model background (from inclusive Z{ee) 
and direct photon production). The dotted line corresponds to the Z{ee) + X 
background alone. Adapted from [25]. 



any new particle via this mechanism is large and depends only weakly on the particle 
mass, in contrast with an exponential dependence characteristic of direct production. 

A relatively large fraction of prompt and energetic photons, electrons, and muons 
expected in the high-multiplicity black hole decays would make it possible to select 
pure samples of black holes, which are also easy to trigger on [25, 26] . At the same time, 
only a small fraction of particles produced in the black hole decays are undetectable 
gravitons and neutrinos, so most of the black hole mass is radiated in the form of 
visible energy, making it easy to detect. 

It has been argued [36] that the fragmentation of quarks and jets emitted in the 
black hole evaporation might be significantly altered by the presence of a dense and 
hot QCD plasma ("chromosphere") around the event horizon. If this argument is 
correct, one would expect much softer hadronic component in the black hole events. 
However, we would like to point out that one would still have a significant number of 
energetic jets due to the decay of weakly interacting W/Z and Higgs bosons, as well 
as tau leptons, emitted in the process of black hole evaporation and penetrating the 
chromosphere before decaying into jetty final states. In any case, tagging of the black 
hole events by the presence of an energetic lepton or a photon and large total energy 
deposited in the detector is a fairly model-independent approach. 

In figure 3 we show the number of black hole events tagged by the presence of an 
energetic electron or photon among the decay products in 100 fb~^ of data collected at 
the LHC, along with estimated backgrounds, as a function of the black hole mass [25]. 
It is clear that very clean and large samples of black holes can be produced at the LHC 
up to the Planck scale of ~ 5 TeV. Note that the black hole discovery potential at the 
LHC is maximized in the e/ 11 + X channels, where background is much smaller than 
that in the 7 + X channel (see figure 3). The reach of a simple coimting experiment 
extends up to Mq « 9 TeV (n = 2-7), for which one would expect to see a handful of 
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black hole events with negligible background. 

The lifetime of a black hole can be estimated using Stefan's law of thermal 
radiation. Since black hole evaporation occurs primarily in three spatial dimensions, 
the canonical 3-dimensional Stefan's law applies, and therefore the power dissipated 
by the Hawking radiation per unit area of the event horizon is p = crT^, where a 
is the Stefan-Boltzmann constant and Th is the Hawking temperature. Since the 
effective evaporation area of a black hole is the area of a 3D-sphere with radius Rs 
and Stefan's constant in natural units is (j = 7r^/60 ^ 1, dropping numeric factors 
of order unity we obtain the following expression for the total power dissipated by a 
black hole: P ~ i?|r^ - Rg^. 

The black hole lifetime r then can be estimated as: r Mbh/P ~ Mbh -R|, and 
using (2) we find: 



Therefore, a typical lifetime of a mini black hole is ~ 10~^^ — 10~^® s. A multi-TeV 
black hole would have a relatively narrow width ~ 100 GeV, i.e. similar to, e.g., a W 
or Z' resonance of a similar mass. This is not surprising, as the strength of gravity 
governing the black hole evaporation rate is similar in the model with large extra 
dimensions to that of electroweak force responsible for the W or Z' decay rates. 

7.3. Accounting for the black hole angular rn,omentum, and grey-body factors 

In the above discussion we used a number of approximations in deriving production 
cross section and decay properties of mini black holes. While reliable accounting for 
more complicated effects related to quantum properties of black holes (spin, quantum 
numbers, etc.) is not possible without intimate knowledge of the underlying theory of 
quantum gravity, a number of authors attempted to estimate some of the above effects 
using simple semiclassical approach. While we don't believe that these estimates are 
any more reliable than the ones obtained in the above simple approximation, we discuss 
these refinements here in some detail. 

The most studied properties of mini black holes beyond our simple picture are 
the effects of its angular momentum and grey-body factors, which have to do with the 
emissivity of particles of various types in the process of black hole evaporation. 

An emissivity of a certain type of particle depends, in general, on its spin, 
S. Indeed, for a small black hole the Schwarzschild radius is comparable with the 
Compton wavelength of the emitted particles. Thus the wave function of a spin-0 
particle would have more overlap with the black hole event horizon than that for 
a spin-1/2 or spin-1 particle. Consequently, one would expect that scalar particles 
are emitted by a black hole more efficiently than spin-1/2 or higher-spin particles. 
This qualitative feature can be parameterized as a grey-body factor, g{S,x), where 
X = E/Th is a dimensionless variable proportional to the energy of the particle. Thus, 
the modified expression (5) for the emitted flux is: 



dx e^ + c 

The grey-body factors can be calculated classically in general relativity. This 
calculation has been extended recently to the case of multidimensional black holes in 
the model with large extra dimensions by a number of authors [37, 38, 39]. 



+ 3 




(8) 



df g{S,x)x^ 



(9) 
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In special cases, the grey-body factors can be calculated analytically, but most 
general calculations that exist to date are performed numerically. The grey-body 
factor for the emissivity of spin- 1/2 particles as a function of their energy is given in 
figure 4 from [38]. As seen from the figure, the effect of grey-body factors is small for 
characteristic black-body radiation energies and n > 2. 
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Figure 4. Grey-body factors for the emissivity of spin-1/2 particles on the brane 
as a function of their energy, E. Vertical axis shows the effective area of the black 
hole horizon for the emission of a spin-1/2 particle expressed in the units of the 
geometrical surface of the event horizon of a three-dimensional black hole, AttR'^. 
Adapted from [38]. 



Recently, grey-body factors for the brane and bulk emission of gravitons have 
been calculated as well [39]. While complete calculations applicable to all graviton 
energies are not available yet, it appears that the fraction of Hawking radiation of bulk 
and brane gravitons is relatively small, except for the cases of very high black-hole 
angular momentum or very large number of extra dimensions (n = 7). 

Since black holes in particle collisions are produced by particles with non- 
vanishing impact parameter, they may carry non-zero angular momentum and thus 
produce a spinning black hole. Given that the characteristic impact parameter is Rs, 
while the relative momentum of the colliding particles is Mbh, one expects a black 
hole produced in such a grazing collision to carry angular momentum 

,'AfBH\^ 

-fl5i\JBH 



Md 

For black holes with masses close to the fundamental Planck scale, typical angular 
momentum is of the order of one. Note that although numerically this is rather small 
angular momentum, it is close to the maximum angular momentum that a black hole 
of such a small mass could have. 

For a rotating black hole, the solution of Einstein's equation is given by Kerr's, 
rather than Schwarzschild's formula and contains explicit dependence on the angular 
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momentum of the black hole. A number of authors have studied properties of Kerr 
black holes in theories with extra dimensions [40]. Both the spin-down effects and 
modification of the grey-body factors due to the angular momentum of the black hole 
have been looked at. In general, for a large number of extra dimensions and for a 
large initial angular momentum of the black hole it is expected that modification 
of the classical black-body radiation picture becomes non-trivial and will have to be 
taken into account for an accurate description of black-hole evaporation. 

Several authors looked at other details of Hawking evaporation process. The 
effects of using the microcanonical ensemble approach, which takes into account that 
the energy of the emitted particles is comparable to the black-hole mass, have been 
discussed in [41] and generally result in the increased black-hole lifetime. The recoil 
effect in the evaporation has been studied [42] as well. 

To summarize, while the recent angular momentum and grey-body factor results 
arc important and encouraging, for the black holes with the masses close to the 
fundamental Planck scale they are likely to be modified in a profound way via unknown 
quantum corrections. Thus detailed studies of the particle content in black- hole 
evaporation probably won't be possible until either their discovery at the LHC or 
a formulation of complete theory of quantum gravity. 

7.4- Other black hole signatures and studies 

While the detection of Hawking radiation from rapidly evaporating TcV-scale black 
holes remains one of the most clear and well-studied signatures at the LHC, several 
authors studied complementary signatures as well. Among those most notable are 
studies of the suppression of dijet production at the invariant masses exceeding the 
fundamental Planck scale [23, 43], anomalous inclusive hadron production due to 
products of black hole decays [44] , and production of black holes in heavy ion collisions 
at the LHC [45]. 

Recently, there have been attempts to estimate the decrease of the black-hole 
production cross section due to inelasticity in parton collisions. The reason for the 
inelasticity is the emission of gravitational waves during the formation of black holes, 
some of which may not be trapped within the event horizon and escape. For discussion 
of this effect, see [46]. 

7.5. Testing Hawking radiation in decays of black holes 

A sensitive test of properties of Hawking radiation can be performed by measuring 
the relationship between the mass of the black hole (reconstructed from the total 
energy of all the visible decay products) and its Hawking temperature (measured 
from the energy spectrum of the emitted particles). One can use the measured 
Mbh vs. Th dependence to determine both the fundamental Planck scale Mjj 
and the dimensionality of space n. This is a multidimensional equivalent of the 
Wien's law of thermal radiation. It is particularly interesting that the dimensionality 
of extra space can be determined in a largely model-independent way via taking 
a logarithm of both parts of (4) for Hawking temperature: logiQ(Tff/lTey) = 
— ^jij- \ogiQ{MBii/^TeV) + const, where the constant does not depend on the black 
hole mass, but only on the Mo and on detailed properties of the bulk space, such as 
relative size of extra dimensions [25]. Therefore, the slope of a straight-line fit to the 
logiQ{TH/lTeV) vs. logiQ{MBu/^TeV) data offers a direct way of determining the 
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dimensionality of space. The reach of this method at the LHC is illustrated in figure 5 
and discussed in detail in [25]. Note that the determination of the dimensionality of 
space by this method is fundamentally different from other ways of determining n, 
e.g. by studying a monojet signature or a virtual graviton exchange processes, also 
predicted in the models with large extra dimensions. The latter always depend on 
the volume of extra space, and therefore cannot provide a direct way of measuring n 
without making assumptions about the relative size of large extra dimensions. The 
former depends only on the area of the event horizon of a black hole, which is not 
sensitive to the size of large extra dimensions. 
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Figure 5. Determination of the dimensionality of space via Wien's displacement 
law at the LHC with 100 ib'^ of data. Adapted from [25]. 



An interesting possibility studied in [27] is production of a precursor of a black 
hole, i.e. a long and jagged highly exited string state, dubbed as a "string ball" 
due to its folding in a ball- like object via a random walk. As shown in [27], there 
are three characteristic string ball production regimes, which depend on the mass of 
the string ball M, the string scale Ms < Mjj, and the string coupling gg < 1. For 
Ms < M < Ms/gs, the production cross section increases oc M^, until it reaches 
saturation at M ~ Ms/gs and stays the same up to the string ball mass ~ Ms/ gl, at 
which point a black hole is formed and the production cross section agrees with that 
from [25]. 

A string ball has properties similar to those of a black hole, except that its 
evaporation temperature, known as Hagedorn temperature [47], is constant: Ts = 
Ms/{2^/2tt). Thus, the correlation between the temperature of the characteristic 
spectrum and the string ball mass may reveal the transition from the Hagedorn to 
Hawking regime, which can be used to estimate Ms and gs- Another possibility is a 
production of higher-dimensional objects, e.g. black p-branes, rather than spherically 
symmetric black holes {p = 0) [48]. For a detailed review see, e.g. [49]. 



7. 6. Simulation of black hole production and decay 

In order to study properties of black holes at colliders, it's important to have tools 
capable of simulating their production and decay. Currently, there are two Monte 
Carlo generators capable of doing this: TRUENOIR [50] and CHARYBDIS [51]. 
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Figure 6. (a) A typical black hole event, generated with CHARYBDIS [51], as 
seen by the ATLAS detector (courtesy of ATLAS- Japan group), (b) A typical 
black hole event, generated with TRUENOIR [50], as seen by the CMS detector 
(courtesy of Albert De Roeck). 



The TRUENOIR generator [50], available since 2001, is a plug- in module for the 
PYTHIA [52] IVIonte Carlo package. It simulates the production and decay of black 
holes as described in [25] , assuming conservation of the individual baryon and lepton 
numbers, as well as the electric and colour charges. It does not include any grey-body 
factors for the black-hole decay and assumes that the evaporation takes place in a 
single, Schwarzschild phase, thus ignoring any spin effects. While the deficiency of 
these assumptions is obvious, the logic behind this simple approach is that given the 
uncertainties related to unknown quantum corrections, an incremental improvement 
from taking into account classically calculated corrections is minor. 

The CHARYBDIS generator [51] is available since 2003 and is interfaced with 
the HERWIG [53] Monte Carlo package. The authors of CHARYBDIS incorporated 
classically calculated grey-body factors and also gave the user flexibility of defining 
when the Schwarzschild stage terminates and a stable Planckian remnant of the black 
hole is formed. 

Both generators have been successfully interfaced with the LHC detector 
simulation packages. Figure 6 shows simulations of a black hole event in the ATLAS 
and CJVIS detectors, using the CHARYBDIS and TRUENOIR JVIonte Carlo packages, 
respectively. A number of dedicated studies [54] of the black-hole decay signatures at 
the LHC have been performed using these JVIonte Carlo generators. 

7. 7. Discovering new particles in black hole decays 

As was mentioned earlier, new particles with the mass ^ 100 GeV would be produced 
in the process of black hole evaporation with a relatively large probability: ^ 1% 
times the number of their quantum degrees of freedom. Consequently, it may be 
advantageous to look for new particles among the decay products of black holes in 
large samples accessible at the LHC and other future colliders. 

As an example, following [55], we demonstrate the discovery potential of a black- 
hole sample to be collected at the LHC for a light Higgs boson. We pick the Higgs 
boson mass of 130 GeV, which is still allowed in low-scale supersymmetry models, 
but makes it quite hard to discover Higgs via direct means either at the Fermilab 
Tevatron [56] or at the LHC. We consider the decay of the Higgs boson into a pair of 
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jets (with the branching fraction of 67%), dominated by the bb final state (57%), with 
an additional 10% contribution from the cc, gg, and hadronic tt final states. 

We model the production and decay of the black holes with the TRUENOIR 
Monte Carlo generator [50]. We used a 1% probability to emit the Higgs particle in 
the black hole decay. We reconstruct final state particles within the acceptance of 
a typical LHC detector and smear their energies with the expected resolutions. We 
select the black hole events by requiring an energetic electron, photon, or a muon, as 
well as total multiplicity of energetic objects (jets, electrons, photons, or muons) of at 
least four. 

The simplest way to look for the Higgs boson in the black hole decays is to use 
the dijet invariant mass spectrum for all possible combinations of jets found among 
the final state products in the above sample. This spectrum is shown in figure 7 for 
Md = 1 TeV and n = 3. The three panes correspond to all jet combinations (with 
the average of approximately four jet combinations per event), combinations with at 
least one 6-tagged jet, and combinations with both jets 6-tagged. (We used typical 
tagging efficiency and mistag probabilities of an LHC detector to simulate 6-tags.) 

The most prominent feature in all three plots is the presence of three peaks with 
the masses around 85, 130, and 175 GeV. The first peak is due to the hadronic decays 
of the W and Z bosons produced in the black hole decay either directly or in top-quark 
decays. (The resolution of a typical LHC detector does not allow to resolve the W 
and Z in the dijet decay mode.) The second peak is due to the h jj decays, and 
the third peak is due to the t — > Wb — > jjb decays, where the top quark is highly 
boosted. In this case, one of the jets from the W decay sometimes overlaps with the 
prompt 6-jet from the top quark decay, and thus the two are reconstructed as a single 
jet; when combined with the second jet from the W decay, this gives a dijet invariant 
mass peak at the top quark mass. The data set shown in figure 7 corresponds to 50K 
black hole events, which, given the 15 nb production cross section for Md = 1 TeV 
and n = 3, is equivalent to the integrated luminosity of 3 pb~^, or less than an hour 
of the LHC operation at the nominal instantaneous luminosity. The significance of 
the Higgs signal shown in figure 7a is 6.7a, even without 6-tagging involved. 

With this method, a 5a discovery of the 130 GeV Higgs boson may be possible 
with £ w 2 pb-i (first day), 100 pb^^ (first week), 1 fb-^ (first month), 10 fb"! 
(first year), and 100 fb~^ (one year at the nominal luminosity) for the fundamental 
Planck scale of 1, 2, 3, 4, and 5 TeV, respectively, even with an incomplete and poorly 
calibrated detector. If the Planck scale is below < 4 TcV, the integrated luminosity 
required is significantly lower than that for the Higgs discovery in direct production. 

While this study was done for a particular value of the Higgs boson mass, the 
dependence of the new approach on the Higgs mass is small. Moreover, this approach 
is applicable to searches for other new particles with the masses ~ 100 GeV, e.g. low- 
scale supersymmetry [57] . Light slepton or top squark searches via this technique may 
be particularly fruitful. Very similar conclusions apply not only to black holes, but to 
intermediate quantum states, such as string balls [27], which have similar production 
cross section and decay modes as black holes. In this case, the relevant mass scale is 
not the Planck scale, but the string scale, which determines the Hagedorn evaporation 
temperature. 

Large sample of black holes accessible at the LHC can be used even to study some 
of the properties of known particles, see, e.g. [58]. 
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Figure 7. Dijet invariant mass observed in black hole decays with a prompt 
Icpton or photon tag in Ri3 pb^-"- of the LHC data, for Mo = 1 TcV and n = 3; 
(a) all jet eombiiiations; (b) jet combinations with at least one of the jets tagged 
as a fe-jet; (c) jet combinations with both jets tagged as 6-jets. The solid line is 
a fit to a sum of three Gaussians and a polynomial background (shown with the 
dashed line). The three peaks correspond to the W/Z bosons, the Higgs boson, 
and the top quark (see text). Prom [55]. 



7.8. Black holes in cosmic rays 

Soon after the original papers [25, 26] on production of black holes at accelerators 
have appeared, it has been suggested that similar black hole production can be 
also observed in the interactions of ultra-high-energy neutrinos with the Earth or 
its atmosphere [59]. For neutrino energies > 10^ GeV, the black hole production 
cross section in uq collisions would exceed their SM interaction rate (see figure 8). 
Several ways of detecting black hole production in neutrino interactions have been 
proposed [59, 60, 61, 62], including large-scale ground-based arrays, space probes, and 
neutrino telescopes as detecting media. 

If the fundamental Planck scale is sufficiently low (1-3 TeV), up to a hundred of 
black hole events could be observed by, e.g. Pierre Auger observatory even before the 
LHC turns on. These estimates are based on the so-called guaranteed, or cosmogenic 
neutrino flux, see, e.g. [63]. In certain cosmological models, this flux could be 
signiflcantly enhanced by additional sources of neutrino emission, e.g. active galactic 
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Figure 8. Cross sections a{uN -> BH) for = Mg™ = 1 TeV and 

n = 1, . . . , 7. (The last four curves are virtually indistinguishable.) The dotted 
curve is for the standard model process i/N — » £X. From [59]. 



nuclei; in this case even larger event count is possible. 

There are two ways to tell the neutrino interaction that results in a black hole 
formation from the standard model processes. The first is based on a particular 
particle content in the black hole events, and would require good particle identification, 
perhaps beyond the capabilities of the existing detectors. The second approach is 
based on the comparison of the event rate for Earth-skimming neutrinos (i.e., those 
that traverse the Earth crust via a short chord, close to the surface) with that for the 
quasi- horizontal neutrinos (i.e., those that do not penetrate the Earth, but traverse 
the atmosphere at a small angle). In the former case, many of the neutrinos would 
be stopped in the Earth due to the large cross section of black hole production. That 
would suppress the rate of the Earth-skimming-neutrino events in a typical ground 
array detector, such Pierre Auger. At the same time, the rate of the quasi-horizontal 
events would increase, as the total cross section, which governs this rate, is dominated 
by black hole production and therefore is higher than in the standard model case. 
By measuring the ratio of the two rates, it is possible to distinguish the standard 
model events from black hole production even with a handful of detected events [61]. 
As model-dependent as they are, current limits on the fundamental Planck scale 
derived [64] from ultra-high-energy cosmic ray data, are already comparable with 
those from sub- millimeter gravity measurements [6] and colliders [8], see figure 9. 

Another interesting observation can be done with the IceCube large-scale neutrino 
telescope at the South Pole, by looking at the zenith angle dependence of the neutrino 
events at various incidental energies. Similar to the previous argument, significant 
reduction of the number of observed events due to the neutrino absorption via black 
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Figure 9. Bounds on the fundamental Planck scale from tests of Newtons law on 
sub-millimeter scales, bounds on supernova cooling and neutron star heating, 
dielectron and diphoton production at the Tevatron, and non-observation of 
black-hole production by cosmic neutrinos. For details and estimates of the 
uncertainties, see [64]. 



hole production in the Earth material surrounding the detector, would occur at smaller 
zenith angles than that in the case of the standard neutrino interactions [62]. In 
addition, particle identification capabilities of the IceCube detector are likely to make 
it possible to detect the black hole events directly by looking at the event shape. 
While it has not been mentioned in the original papers [62], we would like to note 
an additional azimuthal dependence of the event rate for high-energy Earth-skimming 
neutrinos in the IceCube due to the presence of several mountain ridges near the South 
Pole (particularly, the Transantarctic Mountain Ridge). These mountains are not thick 
enough to significantly reduce the fiux of Earth skimming high-energy neutrinos due 
to the standard model interactions, but are sufficiently thick to absorb these neutrinos 
if the black hole creation is allowed. That would constitute a spectacular signature. 

8. Randall-Sundrum black holes 

Mini black holes can be also produced in TeV particle collisions in the Randall- 
Sundrum model [65, 66, 67]. In this case, the warp-factor-suppressed Planck scale, 
~ 1 TeV, plays the role of the fundamental Planck scale in the model with large 
extra dimensions. 

The event horizon of the Randall-Sundrum black holes has a pancake shape with 
the radius in the fifth dimension suppressed compared to the radius Rs on the standard 
model brane by the warp factor e~'^'^^'=. Thus, for Rse^'"^^" ^ nRc^ the black hole 
can be considered "small" and has properties similar to that in the n = 1 (5D) large 



Black holes at future colliders and beyond 



24 



extra dimensions scenario, if the effects of the curvature of the AdS space at the 
standard model branc arc ignored. 

In order to derive properties of the Randall-Sundrum black holes, it is convenient 
to introduce the fundamental 5D Planck scale M, which enters the Lagrangian of 

the Randall-Sundrum model. The relationship between the reduced 4-dimensional 
Planck scale Mpi and M is as follows: Mpi = ^(1 - e-^''''"") « M^/k. Since 
k is 0.01-0.1 xMpi, M = 0.2-0.5 Mpi, i.e. both the 5D and 4D Planck scales are 
of the same order. Since the curvature of the slice of the AdS space is given by 
k'^/M'^ ^ fc^ <§; 1 [66, 67], one can indeed ignore higher-order curvature effects and 
consider Randall-Sundrum black holes as if they were black holes in flat Minkowski 
space. 

The Schwarzschild radius of a black hole of mass Mbh is given by [22, 28, 66] ||: 



Rs ' I 



Taking into account w fcMpi = A^fce^'^'^^s we get: 

Rs « ^,1^. (10) 
VSttA^ y kA„ 

Since the expression under the square root is ~ 10 for a typical range of A'/bh/A^ = 
0(1) and k/Mp\ = 0(0.01), we find that a typical Schwarzschild radius of the Randall- 
Sundrum black hole is Rs ^ 1/A TT 1 TeV , similar to that for the black holes in 
models with large extra dimensions. Indeed, using (2) for the ADD model with n = 1, 
we get: 



i?s (ADD, 5D) 



8Mbh 



which turns into (10) for A^r = Mjj and k = l/8n ^ 0.04. 

Moreover, it is easy to see that such a black hole is still small from the point of 
view of the 5th dimension, as the condition of the black-hole "smallness" mentioned 
above can be expressed as: 

„ nRc knRc kwRc 36 

Rs < = = ~ — -. 

exp(— TrfciJc) ^Mpi exp(-7rA;i?c) A^fc A^fc 

Given that k is between 0.01 and 0.1, the inequality becomes: 

360-3600 
Rs^ — 7 , 

which is clearly satisfied for Rs ^ 1/A^. In fact, one would need to produce a black 
hole with the mass ~ 10^ TeV to exceed this limit. Such energy is achievable neither 
at any foreseen collider nor in fixed-target collisions of ultra-high-energy particles from 
cosmic accelerators. 

II Note that this expression differs from the analogous expression (9) in [65] by a V2 factor; the 
difference stems from the fact that the mass parameter M used in [65] is different from the true 5D 
Planck scale, which enters in the Lagrangian of the model, which we call M in this review. 
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Hawking temperature of the Randall-Sundrum black hole can be found from 
expression (4) for a black hole in models with large extra dimensions by requiring 
n = 1, i.e. 

= 2^' ^''^ 

which, given Rs ^ l/A^r makes it very similar to that for the case of large extra 
dimensions. Consequently, for the preferred range of model parameters both the 
production cross section and the decay properties of a Randall-Sundrum black hole 
are very similar to those in models with large extra dimensions. In fact, both the 
TRUENOIR [50] and CHARYBDIS [51] generators can be used to simulate Randall- 
Sundrum black holes by setting n = l and Md = A„ \/ 8nk « 0.765Mik~'^/^ . 

There has been a suggestion [68] that black holes in Randall-Sundrum model 
and in large extra dimensions can nevertheless be distinguished by the different 
dynamics of an early stage of black-hole evaporation due to the fact that the angular 
momentum of a Randall-Sundrum black hole, unlike that for a black hole in large 
extra dimensions, cannot have any bulk component (due to the existence of a discrete 
Z2 orbifold symmetry). Thus, bulk evaporation for a Randall-Sundrum black hole is 
suppressed compared to that for a black hole of a similar mass in models with large 
extra dimensions. Since it is argued that the bulk component of Hawking radiation of 
gravitons for a black hole in large extra dimensions may be significant during the early 
stages of its evaporation, it is suggested that the black-hole evaporation may result 
in less missing energy in the Randall-Sundrum scenario. Of course, there are other 
ways of distinguishing the two cases, particularly the excitation of narrow graviton 
TeV-scale resonances in Drell-Yan and diboson production in the Randall-Sundrum 
case [13], versus an overall enhancement of the high end of the Drell-Yan and diboson 
mass spectra in the case of large extra dimensions. 

Recently, there have been studies of modification of black-hole properties due 
to Gauss-Bonnet or Lovelock terms added to the Einstcin-Hilbert action. These 
modifications affect properties of black holes in models with either large or warped 
extra dimcnsion(s). Additional terms could naturally introduce a minimum threshold 
on the black-hole mass. For detailed studies of modifications related to these higher- 
order terms, see [66, 67]. 



9. Black holes at CLIC 



A discovery of black holes at the LHC or in ultra-high-energy cosmic ray interactions 
would open a whole new avenue of studies of quantum gravity. If this is the case, it 
would be very important to accumulate large samples of black holes in well-controlled 
conditions, with various initial quantum numbers, including black holes with the 
quantum numbers of vacuum. 

Since most of black holes produced at the LHC would carry electric charge 
and baryon number, as well as broad range of masses, it is important to have a 
complementary way of producing completely neutral black holes of fixed mass. An 
e+e~ machine of sufficient energy would open such possibility^. 

Given the existing limits on the TeV-scale gravity, it is clear that such a 
machine would have to have multi-TeV energy, i.e. beyond the reach of the proposed 

^ A multi-TeV muon collider would offer similar capabilities, although technologically such a machine 
is more far-fetched. 
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Figure 10. Black hole properties at a 5 TcV CLIC e+e^ collider. Plots a)-c) 
correspond to the production cross-section, Hawking temperature, and average 
decay multiplicity for a fixed-mass 5 TeV black hole. Plot d) shows differential 
cross section of black hole production for n = 4, as a function of its mass. 



International Linear Collider. The only known technology that would allow to bnild 
such a machine is two-beam acceleration, which is the basis of a proposed CLIC 
accelerator with the energy reach of 3-5 TcV [69]. For completeness, we include 
properties of black holes produced at a 5 TeV CLIC collider in this review. They are 
shown in figure 10. 

There are several advantages of producing black holes in e+e~ collisions. First 
of all, unlike that at hadron colliders, nearly all black holes will be produced with 
the same mass, equal to the total energy of the machine. The spread due to the 
"beamstrahlung" is expected to be small, as shown in figure lOd. Second, the black 
holes are produced at rest and have the vacuum quantum numbers. That would allow 
to study the details of evaporation process. Third, since the energy of the e~^e~ 
collider can be readily tuned, one could produce black holes close to the quantum 
production threshold {~ Mo ) and study the effects of quantum gravity near the 
threshold. Finally, using polarized beams one can change the angular momentum of 
produced black holes and study how it affects their evaporation. The ability to study 
all of these fine effects would make CLIC a very attractive complement to the LHC 
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for studies of black holes. 

10. Reentering black holes 

An interesting topic in black hole phenomenology, which has not been studied in much 
detail yet, is the possibility that a black hole, once produced, moves away into the bulk 
space. Normally it does not happen as the black holes produced in collisions at the 
LHC or in cosmic ray interactions are likely to have charge, colour, or lepton/baryon 
number hair that would keep them on the brane. However, a possibility of that kind 
is allowed in the case when the strength of gravity in the bulk and on the brane is 
very different. This is the case, e.g. in the scenario with large extra dimensions with 
an additional brane term [70] , or in the case of infinite- volume extra dimensions [71] . 

In these models, a particle produced in a subplanckian collision, e.g. a graviton, 
coiild move away in the bulk, where it bcx'onies a })lac;k hole due to much lower effective 
Planck scale in extra dimensions. Since the Planck scale in the bulk is very low, e.g. 
~ 0.01 eV in the infinite- volume scenario [71], the newly- formed black hole is very 
cold and therefore essentially stable. Furthermore, it generally does not move far 
away from the brane due to gravitational attraction to it, and can further accrete 
mass from relic energy density in the bulk and from other particles produced in the 
subsequent collisions. Once the mass of the black hole reaches the mass of the order of 
the apparent Planck scale, Mpi ~ 10^^ GeV, the event horizon of the bulk black hole 
grows so large that it touches the brane, and the black hole immediately evaporates 
on the brane into ^ 10 particles with the energy ^ lO'^^ GeV each. (The energy 
released in such an event is similar to that in an explosion of a large, few hundred 
pound conventional bomb!) If such black holes are copiously produced by a remote 
cosmic accelerator of a reasonable energy, they could act as a source of the highest- 
energy cosmic rays that are emitted in the process of decay and deceleration of the 
super-energetic black hole remnants. 

Even if the mass of the black hole in the bulk is small, it has certain probability 
to reenter our brane. In this case, since the event horizon cannot be destroyed, once 
it has been formed, such a subplanckian object would likely to act as a black hole on 
the brane and evaporate similarly to a transplanckian black hole discussed above. 

11. Conclusions 

To conclude, black hole production at the LHC and in cosmic rays may be one of 
the early signatures of TeV-scale quantum gravity in the models with large or warped 
extra dimensions. It has three advantages: 

(i) Large cross section: because no small dimcnsionlcss coupling constants, analogous 
to a, suppress the production of black holes; this leads to enormous rates. 

(ii) Hard, prompt, charged leptons and photons: because thermal decays are flavour- 
blind; this signature has practically vanishing standard model background. 

(iii) Little missing energy: because most of the black hole evaporation products are 
detectable; this facilitates the determination of the mass and the temperature of 
the black hole, and may lead to a test of Hawking radiation. 

Large samples of black holes accessible by the LHC and the next generation of 
colliders would allow for precision determination of the parameters of the bulk space 
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and may even result in the discovery of new particles in the black hole evaporation. 
Limited samples of black hole events may be observed in ultra-high-energy cosmic ray 
experiments, even before the LHC turns on. 

If low-scale gravity is realized in nature, the production and detailed studies 
of black holes in the lab are just few years away. That would mark an exciting 
transition for astroparticle physics: its true unification with cosmology — the "Grand 
Unification" to strive for. 
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